Amplitude death in complex networks induced by environment 
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We present a mechanism for amplitude death in coupled nonlinear dynamical systems on a complex 
network having interactions with a common environment-like external system. We develop a general 
stability analysis that is valid for any network topology and obtain the threshold values of coupling 
constants for the onset of amplitude death. An important outcome of our study is a universal 
relation between the critical coupling strength and the largest non-zero eigenvalue of the coupling 
matrix. Our results are fully supported by the detailed numerical analysis for different network 
topologies. 
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1. INTRODUCTION 

The dynamics of most real world systems is very com- 
plex and can be analyzed by considering them as many 
sub-systems or units interacting with each other. Such 
systems then can be modeled by complex networks of di- 
verse topologies chosen to suit their collective behavior. 
It is now well established that the interaction among the 
sub-systems can lead to emergent phenomena like syn- 
chronization and amplitude or oscillator death [D-Q- In 
such cases, the global dynamics depends on the inter- 
play between the network structure and nodal dynamics. 
However, the interactions of complex networks with en- 
vironment, such as an external agency or medium, and 
the consequent emergent dynamics are not studied yet. 
In the case of quantum systems, the role of environment 
in causing decoherence, relaxation, and dissipation is well 
studied [J-Q . In the context of biological systems, the en- 
vironment can play a constructive role as the mechanism 
for triggering or signaling coordinated rhythms [8Hl0|. 
Hence it is important and relevant to study the effect of 
environment on the dynamics of coupled systems. 

Among the emergent phenomena in coupled systems, 
the suppression of dynamics or amplitude death is of- 
ten a useful control mechanism for stabilizing systems to 
steady states. It is an important self organized behavior 
that can play crucial roles in regulatin g, s witching, and 
controlling physica l fill Il2| . chemical [13l - [l6j . and bio- 
logical systems In this context, amplitude death 
(AD) refers to the phenomenon where the coupled or in- 
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teracting sub-systems settle to a steady state in which 
dynamics is quenched. 

In the context of regularly coupled systems, it has been 
shown that amplitude death can be induced by different 
mechanisms, such as parameter mismatch [15l . |22l \2c 
time-delay coupling [2J-[28(> conjugate coupling [2£ 
attractive and repulsive couplings 

coupling [32]. In addition to the study of amplitude 
death in regularly coupled systems, there has been a 
few recent studies on amplitude death in complex net- 
works. In the specific context of networks, amplitude 
death has been studied for parameter mismatch or detun- 
ing of frequencies in an ensemble of limit cycle oscillators 
with mean field coupling [33| , array of limit cycle oscil- 
lators 34; 36j .small- world networks 21], and scale- free 
networks [33]. So also, time-delay in coupling is found to 
induce amplitude death in networks of limit-cycle oscil- 
lators [13,[3i-[i3 and chaotic systems (Mill- Recently, 
it was shown that it is possible to target amplitude death 
in a network of nonlinear oscillators by a proper choice 
of nonlinear coupling [43^ . 

In this paper, we present an interesting phenomenon 
in which the collective dynamics of coupled systems is 
quenched due to an interaction with an environment or 
external agency. For this, we model the average effect 
of the environment by an over-damped oscillator which 
is kept alive with feedback from the subunits. We find 
that while the coupling among the units can give rise to 
a synchronizing tendency, the coupling through the envi- 
ronment has a tendency to drive the systems to a state 
where the sum of the variables is small. The combined 
effect of these two tendencies is to lead the coupled sys- 
tems to the state of amplitude death. Our method has 
the advantage that it involves a single damped dynamical 
system coupled to all nodes equally and hence the design 
procedure is simple and easy to implement. It is found to 
be effective in complex networks of different topologies. 
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We have found this mechanism to be quite general and 
effective in inducing amplitude death in two systems cou- 
pled by different types of coupling and of different intrin- 
sic dynamics fiij . The present study generalizes the pre- 
vious results in two directions. First, we consider a com- 
plex network of TV systems and develop the stability anal- 
ysis following the approach given in Ref. [45[ . The stabil- 
ity conditions are obtained for the general case. Second, 
we consider different network structures. Our results are 
supplemented by detailed numerical analysis where in- 
dices of amplitude death are computed directly from sim- 
ulations. For numerical simulations, we use Rossler as a 
standard system in its chaotic region. However, we have 
tried this method for Landau-Stuart oscillators and the 
Hindmarsh-Rose model of neurons and found it to be ef- 
fective in causing amplitude death. An important finding 
from our study is that, the critical strength of coupling 
needed for amplitude death has a universal relation with 
the largest non-zero eigenvalue of the coupling matrix 
which is tested for many symmetric networks like chain, 
ring, tree, lattice, all-to-all, star, and random topologies. 



2. AMPLITUDE DEATH VIA DIRECT AND 
INDIRECT COUPLING 

We consider the dynamics of TV systems Xi, i = 
1 , 2, . . . , TV, in a network, coupled with two types of cou- 
plings, namely, a direct diffusive coupling and an indirect 
coupling through an environment as an extension of the 
model given in Ref. [44| which gave a general model for 
amplitude death in two coupled systems. The dynamics 
of such a model is given by 



ii = f(xi) + 2J /3Gije d Xj + e e jw 
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where i, j = 1, 2, ...TV. Here, Xi represents m-dimensional 
nonlinear oscillators whose intrinsic dynamics is given by 
f(xi). G is the coupling matrix of dimension TV x TV. 
We choose the elements of G such that, the row-sum, 
Y/V Gij = 0, for every j. This ensures that the largest 
eigenvalue of the coupling matrix /ii, is zero, ft is a ma- 
trix (mxm) with elements and 1 and defines the compo- 
nents of Xi which take part in the coupling. For simplic- 
ity, we take /3 to be diagonal, (i = diag((3\, /3 2 , . . . , (3 m ) 
and in numerical simulations, only one component, f3\ is 
assumed to be non-zero. The environment is considered 
to be a one-dimensional over-damped oscillator w, with 
damping parameter k. It is clear that without feedback 
from the systems, the environment can not remain dy- 
namic and will rapidly settle to a steady state. However, 
the feedback from all the systems keeps it active. All the 
systems, in turn, get feedback from the environment. 7 
is a column matrix (m x 1), with elements or 1, and 
it decides the components of Xi that get feedback from 
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FIG. 1: Time series of the first variables xn of 10 cou- 
pled Rossler systems in an all-to-all coupled network [Eq. |[2)] 
showing amplitude death for (e<j,e e ) = (0.05,0.8). Here, the 
Rossler parameters are a — b — 0.1, c = 18. The damping 
parameter of the environment is taken to be k = 1. 



the environment. j T is the transpose of 7 and decides 
the components of Xi which gives feedback to the envi- 
ronment. The strength of this feedback coupling between 
the systems and the environment is given by e e . 

We illustrate our scheme using a network of cou- 
pled chaotic Rossler systems represented by the following 
equations: 



Xil 

±i2 
Xi3 

w 



= -Xil - Xi3 + £<i 2J Gij x jl + e e 



xn + ax l2 , 

b + x i3 (xn - c), 



TV 



— Xil. 



(2) 



Here, we choose G to be an all-to-all connected network 



of 10 nodes, that is, G^ = 1, if j ^ i and Gy 



-9. We 



find that amplitude death is possible for suitable values 
of coupling strengths. The time series for the amplitude 
death state is shown in Fig. Q] We note that this state 
corresponds to the stable fixed point of Eq. , given by 



Xil = 


r - Vc 2 - 


x i2 = 


-x*Ja, 


x i3 — 


-b/(x*i - < 


w 


— ^eXu/n. 



Aabn/(K - e 2 e a)}/2, 



(3) 



As noted from Eq. ([2]), our model has two types of cou- 
pling. The first is direct diffusive coupling which tends to 
synchronize the systems Xi = Xi = x% = . . . = xjy. We 
find that, the coupling via environment has a tendency 
to decrease the sum, ■ Xi . When both these tendencies 
work together the systems converge to a fixed point. This 
is explicitly seen in the context of two systems coupled 
through the environment |44j|, where the environmental 
coupling reduces the sum to a small value, corresponding 
to antiphase synchronization. 



3 



3. LINEAR STABILITY ANALYSIS 

The stability of the steady state of the network of cou- 
pled systems given in Eq. ([1} can be analysed by writing 
the variational equations formed by linearising Eq. (JlJ as 

j 

z = -K2-|7 T I>, (4) 

i 

where £j and z are small deviations from the respective 
values of a^and y, and /' is the m x m Jacobian matrix. 
Let us introduce the m x N state function (4f| 

S = (&,$,,..., (5) 

Then, Eq. ((4]) for the synchronized state (x\ = x-i = 
. . . = xjy) can be written as 

3 = f'E + f3e d ZG T + e e zT, (6) 
z = -kz- ^7 T X]^' ( 7 ) 

i 

where G T is the transpose of the coupling matrix, and T 
is an m x N matrix, r = (7, 7, . . . , 7) . 
Let e k be an eigenvector of G T such that 

G T e k = h i k e k , (8) 

where [i k is an eigenvalue of G T . Right-multiplying both 
sides of Eq. ^ with e k , we get 

Ee k = f'Sek + ^ k /3e d Ee k + e e zTe k . (9) 

Let 

$ fc = Ee k . (10) 
Then, Eq. ^ can be written as 

= + Mfc/3eA + e e 2re fc . (11) 

We note that e\ = (1,1,..., 1) T is the synchronization 
manifold and $1 = Sei = ^. Since one could write T 
as the product T = jej, Eqs. (TTTj) and ([7]) can be written 

as 

= /'$k + + e e 27ef e fc , (12) 

z = -K«-^ 7 T *i. (13) 

First, we consider the case where G is taken to be a 
symmetric matrix. In this case, the remaining eigenvec- 
tors span an (N — l)-dimcnsional subspacc orthogonal 
to the eigenvector ei . Consequently, this subspace is 
orthogonal to the synchronization manifold. For k = 1, 
Eq. ([12]) becomes 

*i = /'*! + ^zNy. (14) 



Since e, are orthogonal, Te k = 7&fefc = for k ^ 1. 
Therefore, For fc ^ 1, Eq. (fT2"j) reduces to 

*fc = /'*fc+A*fc^*fc- (15) 

We note that Eqs. (jT4j) and (|L3|) are coupled while 
Eq. (fT5|) is independent of the other two. Moreover, 
Eq. (|15[) is equivalent to the master stability equation 
introduced by Pecora and Carrol in Ref . [4q] . There- 
fore, the stability function for any given system will be 
obtained as a function of e d ^ k in the same way. This 
therefore ensures the stability of the synchronized state 
Xi = X2 = X3 = . . . = xn- As noted in the previous 
section, for the amplitude death state to be stable, we 
need one more condition to be satisfied. That is, the 
synchronized state should be a fixed point. For this, the 
eigenvalues of the Jacobian corresponding to the coupled 
system given in Eqs. (fT4"]) and (| L3[) should be negative. 

So far, we have discussed the case where the coupling 
matrix G is symmetric. The same analysis can be ex- 
tended to the asymmetric case as well. In this case, the 
eigenvectors of G are, in general, not orthogonal to the 
synchronization manifold. 

Let e k , for k 5^ 1, be split to two components, one 
parallel and the other perpendicular to e\. That is, 

ek = e k + q\\ei, (16) 

where, is orthogonal to e\. Substituting e k from 
Eq. (jTHJ) in Eq. (|11[) . we find that the dynamics of $1 
and z will again be the same as given in Eqs. (1141) and 
(fT3| . The dynamics of $ k , k > 1, are given by 

3»fe = f'$k + Vkfad^k + e e zN<yq\\ , (17) 

since Fe^ = jef = 0. In principle, the coupled equa- 
tions, Eqs. ([Lip. (fT7|) and (fl~3~|) can be considered as a 
master stability equation in this case also. However, in 
this case it is not practically useful since the master sta- 
bility function will be a function of four parameters e d , 
a, b, and e e , with fx k = a + ib. 

To continue the analysis of the stability of the am- 
plitude death states from Eqs. (JH) , dD i an( l we 
assume that the time average values of /' are approxi- 
mately the same and can be replaced by an effective con- 
stant value a. In this approximation we treat £j's to be 
scalars. This approximation simplifies the problem such 
that only the relevant features remain and is expected to 
give features near the transition. This type of approxi- 
mation was used in Refs. [3, H3, Hi| and it is found to 
describe overall features of the phase diagram reasonably 
well. 

Thus, Eq. (fT3|) becomes 

= a<P k + fi k Pe d <$> k , (18) 
and the corresponding Lyapunov exponent is given by 
Ai = a + (i 2 ed, (19) 
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where /i 2 is the largest Hk for k 7^ 1. 

The Jacobian corresponding to the coupled Eqs. (JTJ 
and fill) is 



J = 



-e e /N — k 



and the eigenvalues are 

(ch-k)± V^K - a) 2 - 4e2 
A 2,3 — o • 



(20) 



For the stability of the amplitude death state, the 
real parts of the eigenvalues should be negative. Thus 
Eq. (jT9]) gives the condition 



a + u 2 e d < 0, 



(21) 



while from Eq. ([2TJ]) we get the following conditions. 
(1) If (k — a) 2 < 4(e 2 — an), A2,3 are complex and the 
condition of stability is 



k > a. 



(22) 



(2) If (k — Hied — a) 2 > 4(e 2 — an), A2,3 are real and the 
stability condition becomes 



k > (a) and e 2 > (are). 



(23) 



Thus, if Eqs. ([2T| and ([22]) or ([23]) are simultaneously 
satisfied, the oscillations can not occur and the systems 
stabilize to a steady state of amplitude death. 

For a given k, a, and \x 2 , the transition to amplitude 
death occurs at critical coupling strengths ea c and e ec 
independent of each other. That is 



and 



€dc = const 



const 



(24) 



(25) 



For different network configurations /i 2 is different and 
the transition occurs at the critical coupling strength 



—a 
M2 ' 



(26) 



In the case where G is asymmetric, using the approx- 
imation /' ~ a as explained above, we can write the 
Jacobian corresponding to Eqs. (fbf]) . ([TT]) . and (TJ3]) as 

a + H2/3td e e zN^fq\\ e e N^q\\ 
J = I a + /i 2 £d e e N 
-e e /N — k 

and the eigenvalues are the same as given in Eq. (|20p . 
Thus, we get the same stability relations as in Eqs. (|2"TT) . 
and 



4. NUMERICAL ANALYSIS 

In this section, we apply our scheme to different net- 
work topologies. First, we apply the scheme of coupling 
introduced in Eq. (fTJ) to the case of regular networks of 
coupled chaotic Rdssler systems. Here, we take the cou- 
pling to be of a diffusive type [Eq. @] . The occurrence of 
amplitude death in the case of a regular all-to-all coupled 
network is illustrated in Fig. [TJ 

To characterize the state of amplitude death, we use an 
index A introduced in the earlier paper 44j . It is defined 
as the difference between the global maximum and global 
minimum values of the time series of the system over a 
sufficiently long interval. The case where A = repre- 
sents the state of amplitude death, while A 7^ indicates 
oscillatory dynamics. The parameter value at which A 
becomes ~ is thus identified as the threshold for onset 
of stability of amplitude death states. 

For a given network topology, the threshold value of 
coupling strengths for the onset of amplitude death is 
given by Eqs. (|2~4"]) and (|2"5"]) . This is verified for the case of 
an all-to-all coupled network of Rdssler systems by direct 
numerical simulations. Using the index A the region of 
amplitude death states are identified in the parameter 
plane of coupling strengths, e e -ed, and is shown in Fig. [5] 
The transition curves from the stability analysis given in 
Eqs. (f2~4"]) and ([2"5]) are also plotted. We see that the 
agreement is good. 

However, we note that exact agreement with stability 
theory as shown in Fig. [2] is seen only for Rdssler type 
nodal dynamics. For a network of Landau-Stuart oscil- 
lators, given by the following equations, 



xn = (1 - x l ix - xf 2 )xn - ujx i2 + £ d ^ Gi i x J 1 + CeW ' 

3 

x l2 = (1 - x 2 n - x 2 i2 )xi 2 + wxa, 

|E-- ( 27 ) 



W = —KW 



similar analysis shows some deviations between theoreti- 
cal and numerical transition curves (Fig. [3]). The reason 
for this is the following. Towards the end of the stabil- 
ity analysis, we have used an approximation of constant 
Jacobian /', which masks the system-specific details of 
the transition, but gives the overall features of the phase 
diagram. Hence the conditions (24) and (25) are approx- 
imate, and one must investigate in a specific case to see 
any departures from them. 

For both Rossler and Landau-Stuart in the amplitude 
death state, the Jacobian /' depends on e e , but not on 
ed- Hence, the condition (24), i.e., e ec = const, obtained 
from Eq. (23) which is derived from Eqs. (13) and (14), 
is independent of td for both Rossler and Landau-Stuart 
as can be seen from Figs. 2 and 3. 

The other condition (25), i.e., Cdc = const, is obtained 
from Eq. (21) which is derived from Eq. (15). Since Eq. 
(15) depends both on and indirectly on e e through the 
Jacobian /', e^ c will now depend on e e for both Rossler 
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FIG. 2: Transition from the region of oscillations (I) to 
the region of amplitude death (II) is shown in the parameter 
plane e e -£d for the coupled Rossler systems. Numerical sim- 
ulations are done on a symmetric, all-to-all coupled network 
of 10 nodes. The points mark the parameter values (e ec ,edc) 
at which the transition to the amplitude death occurs. Solid 
triangles show the transition to amplitude death as is in- 
creased for a constant e e . The horizontal line formed by these 
triangles confirms the stability condition Eq. (|24|l . Similarly, 
circles correspond to transition to the amplitude death state 
as e e is increased for a constant e<j and confirm the stability 
condition of Eq. ([25j . 



and Landau-Stuart. In the case of Rossler networks, the 
Jacobian has a simple structure and the dependence of 
edc on e e is weak giving almost a straight line as in Fig. 2. 
In the case of Landau-Stuart networks, the dependence 
of €dc on e e is a polynomial relation explaining the curve 
obtained for numerical simulations (Fig. [3]) . 

In the context of Landau-Stuart oscillators, there is an 
additional complexity due to bistability with oscillations 
and amplitude death co-existing with different basins. 
This bistability has been reported earlier in the case of 
amplitude death in Landau- Stuart oscillators [2!| |44| . 
Such a bistability does not exist for Rossler systems. In 
Fig. |3l we have used the same set of initial conditions for 
any pair of (e e ,£d) values. For a different set of initial 
conditions, the critical coupling curve in Fig. [3] can shift 
slightly though the general features will remain the same. 

We also verify numerically the criteria for transition 
to amplitude death given in Eq. (|2U)) . n Fig. [H the line 
corresponds to the stability conditions of Eq. (p?5|) and 
the points are square of the critical coupling strength e e , 
as k is varied. As we can see from Fig. 21 the agreement 
is good for larger values of n. However, for small values 
of k, the points deviate from straight line behavior. The 
reason is clear from Eq. (f22j) which gives the lower limit 
on k. 

The nature of the transitions to the state of amplitude 
death is further characterized by fixing one of the param- 



FIG. 3: Transition from the region of oscillations (I) to the 
region of amplitude death (II) is shown in the parameter plane 
£e-£d for 10 coupled Landau systems on a symmetric, all-to-all 
coupled network [Eq. (|27[l ]. The points mark the parameter 
values (e e c,e<ic) at which the transition to the amplitude death 
occurs. Here, the intrinsic parameter of the systems and the 
damping parameter of the environment are chosen to be cu — 
2 and k = 1 respectively. We use the same set of initial 
conditions for any pair of (e e ,ed) values. 
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FIG. 4: Transition from the region of oscillations (I) to the 
region of amplitude death (II) is shown in the parameter plane 
k — £g for the coupled Rossler systems given in Eq. ([2]). Here, 
an all-to-all coupled network of 10 nodes is used. The points 
mark the square of the coupling strength at which the transi- 
tion to the amplitude death state occurs for each value of n. 
Solid curve is a linear fit corresponding to the stability condi- 
tion Eq. ()23[) . The deviation from straight line behaviour for 
small values of k is discussed in the text. The solid curve is a 
linear fit corresponding to the stability condition of Eq. (23). 
The deviation from straight line behavior for small values of 
k is discussed in the text. 
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FIG. 5: (a) The oscillatory part of incoherent energy, E, as 
a function of e e for a fixed value of ed = 0.4 for an all-to-all 
coupled network of 10 Rossler systems, (b) The index < A > 
as a function of e e for a fixed value of = 0.4 for an all-to-all 
coupled network of 10 Rossler systems. As e e is increased, 
we observe a continuous transition to the state of amplitude 
death. 



eters e e or and increasing the other. To characterize 
the transition, we use the oscillatory part of incoherent 
energy, E defined in Ref. j34[. Since the fixed point ob- 
tained in this case is not the origin, the oscillatory part 
of incoherent energy is defined after shifting the origin to 
the fixed point as 

E < £f=l((Sjl - A? + {Xj2 ~ A? + (X j3 Xjf) > 

< Ef =1 ((^i - *i) 2 + K° 2 - 4) 2 + (4» - 4) 2 ) > ' 

(28) 

where, x® 2 , x%) represent the variables (xj±, Xj2, 

Xjs) in the uncoupled case (e e = ed = 0), (x*, arfj, £3) 
represent the fixed point of the coupled system (Eq. [3]) , 
and < . > denotes average over time. In Fig. [Ha), we 
plot E of the coupled system given in Eq. [2] for increas- 
ing e e for a chosen value of ed- Here, the transition from 
oscillatory state to amplitude death state is continuous 
such that, as the coupling strength is increased, E grad- 
ually decreases to zero. However, the calculation of E, 
using Eq. [28] is useful only when the fixed points of the 
coupled system can be calculated analytically. Alterna- 
tively, we can use the index < A >, used to identify states 
of amplitude death in Figs. 12131 and 14} which does not re- 
quire the knowledge of the fixed point. This is shown 
in Fig. EJb). Numerically, we also observe that, at each 
node the sub-systems undergo a reverse period-doubling 
bifurcation to limit cycle before undergoing a transition 
to the amplitude death state, similar to the case of two 
coupled Rossler systems reported earlier [44|. Similar 
transition is observed for the case where e e is kept fixed 
and ed is increased. 

So far, we have presented the results from numerical 
simulations of Eq. ([2]) in an all-to-all coupled network. 
Similar results arc observed for other network configu- 
rations such as chain, ring, tree, lattice, star and ran- 
dom topologies. From numerical simulations of Rossler 
systems coupled in different network topologies, we see 
that, the critical strength of coupling via environment 
e ec is independent of the network topology. On the other 
hand, the critical strength of direct coupling for ampli- 
tude death, e<j c varies with the largest non-zero eigenvalue 
of the coupling matrix G, as given in Eq. (|26[) . To verify 
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FIG. 6: Critical strength of direct coupling edc for amplitude 
death as a function of the largest non-zero eigenvalue of the 
coupling matrix, /12. Here points correspond to values ob- 
tained from numerical simulations, while the line correspond 
to the stability condition in Eq. (|26[) . Open circles repre- 
sent symmetric networks of different topologies such as chain, 
ring, all-to-all-coupled, tree, lattice, star and random. Simi- 
larly, filled triangles represent asymmetric networks of differ- 
ent topologies such as chain, ring, tree, star and random. For 
the asymmetric networks, /i2 is, in general, complex. Hence, 
the real part of fi2 is plotted here. The parameters of the 
Rossler system are the same as that used in Fig. [T] The pa- 
rameters in the coupling terms are e e = 0.8 and k = 1. 

this, we consider symmetric and asymmetric matrices of 
different topologies and sizes. With each network con- 
sidered, the largest non-zero eigenvalue, H2 of the corre- 
sponding coupling matrix G is calculated. The critical 
value of coupling, Cd c is obtained from numerical simula- 
tions of Eq. ([2]) and is plotted against the corresponding 
fi2 in Fig.[S] A universal relation between the critical cou- 
pling strength and largest non-zero eigenvalue of the cou- 
pling matrix, as given by Eq. (|26[) is clearly seen. Similar 
behaviour is seen in the case of coupled Landau-Stuart 
systems on a network, as shown in Fig. [7] A similar in- 
sensitivity of the transition to amplitude death, to the 
network structure is reported in the case of time-delay 
coupled Rossler systems in Ref. [42| where, the small- 
est eigenvalue of the adjacency matrix of the network is 
found to determine the size of the death island. 



5. DISCUSSION 

We report the amplitude death in complex network of 
nonlinear oscillators caused by interactions with a com- 
mon environment. Our method involves a damped en- 
vironment modeled by a single variable coupled to all 
nodes equally. We develop a stability analysis to obtain 
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FIG. 7: Critical strength of direct coupling e<j c for ampli- 
tude death as a function of the largest non-zero eigenvalue of 
the coupling matrix, /12 for Landau-Stuart oscillators. Here 
points correspond to values obtained from numerical simu- 
lations, while the line correspond to the stability condition 
in Eq. |126[| . Open circles represent symmetric all-to-all cou- 
pled networks of different sizes and filled squares represent 
symmetric networks of chain topology of different sizes. The 
parameters in the coupling terms are e e = 4.0 and k = 1. 



the criteria for the onset of amplitude death. The transi- 
tion curves obtained from the stability analysis matches 
well with those obtained from direct numerical simula- 
tions. Moreover, the method introduced here is found 
to work for different network topologies. In the context 
of the two specific nodal dynamics, Rossler and Landau- 
Stuart studied here, wc find that, there exists a universal 
relation which is independent of network topology, be- 
tween the largest eigenvalue of the coupling matrix and 
the critical value of coupling. All the points correspond- 
ing to transition to amplitude death state for different 
network topologies fall on the same curve. This is as 
expected from the stability analysis developed. 

The dynamical mechanism that induces amplitude 
death itself is very interesting, where the environment 
modulates the dynamics in a self-organized way. Since 
amplitude death is brought about by a common variable 
coupled equally to all nodes, the design procedure is sim- 
ple and easy to implement in cases where targeting of 
complex systems to steady state behaviour is desirable. 
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